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Abstract 

The scale dependence of the evolution of photoproduction cross sections with the 
photon-proton centre of mass energy W is studied using low Q 2 < 0.01 GeV 2 e + p 
interactions collected by the HI experiment at HERA. The value of the largest 
transverse momentum of a charged particle in the photon fragmentation region is 
used to define the hard scale. The slope of the W dependence of the cross section 
is observed to increase steeply with increasing transverse momentum. The result is 
compared to measurements of the Q 2 evolution of the W dependence of the virtual 
photon-proton cross section. Interpretations in terms of QCD and in terms of Regge 
phenomenology are discussed. 
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1 Introduction 



At high energies measurements of the energy dependence of real photon-proton cross 
sections are consistently described by the power law a l °^ oc (W 2 ) x , where W is the jp 
centre of mass energy. The power is A ~ 0.08 |TJ, 0] and is found to be universal for all 
photon or hadron interactions with protons ||. The virtual photon-proton cross section 
however, which is at large W related to the structure function F 2 via cr'°* p ~ 4ir 2 a/Q 2 ■ 
F2(W 2 , Q 2 ), is found to rise fast with increasing W 2 for Q 2 values larger than a few GeV 2 
0, ||. Here Q 2 is the virtuality of the photon and a the fine structure constant. At large 
W the relation W 2 = Q 2 /x holds, with x the deep inelastic scattering (DIS) variable 
Bjorken-x. In the HERA regime fits of the form <r* < ? oc (W 2 ) x lead to values of A which 
rise from 0.2 to 0.4 in the Q 2 range from 1.5 to 10 3 GeV 2 j§. 

In Regge theory the energy dependence of the total cross sections is specified by the 
intercept, which is simply 1 + A, of the pomeron, the leading exchanged trajectory. It is 
argued (see e.g. f| f?j) that the change of slope parameter A with increasing Q 2 is asso- 
ciated with a transition from a non-perturbative 'soft' to a perturbative 'hard' regime. 
Indeed, in the domain where Q 2 is larger than a few GeV 2 and perturbative QCD can 
be safely applied, it predicts an increase of A with Q 2 . In the domain of 'soft' interac- 
tion phenomenology the change of the slope parameter can be interpreted in Reggeon 
Field Theory as a reduction of screening corrections with increasing Q 2 M. Presently 
Reggeon Field Theory and perturbative QCD are two complementary approaches which 
successfully describe physics processes in different regimes, but the transition region be- 
tween these soft and hard regimes is poorly understood. HERA provides an ideal testing 
ground to study the transition from soft to hard physics by comparing real and virtual 
photon-proton collisions in the same experiment. Comparisons of the energy dependence 
of the cross section or of the final states properties (e.g. multiplicities, energy flows, in- 
clusive particle distributions) allow more insight to be gained into the dynamics of strong 
interactions and can be helpful in the development of a common underlying picture in 
terms of an effective field theory (see e.g. ||). 

In a previous publication energy flow in the photon fragmentation region was stud- 
ied [JTOJ . In this paper we investigate the energy behaviour of the cross section for photo- 
production events containing charged particles with high transverse momentum, p t , and 
compare it with the energy behaviour of the total inclusive 7*p cross section. The pres- 
ence of a high p t charged particle allows a scale to be identified which can be varied in 
a similar way as the scale in virtual photon interactions, and the effect on the energy 
dependence of the cross section can be studied. A continuous coverage of both the soft 
and hard scattering domains becomes possible. To emulate the kinematical configuration 
of a 7*p collision as closely as possible the high p t charged track is required to be in the 
photon fragmentation region, which corresponds at the same time to the largest accep- 
tance region of the HI tracking detectors. The idea of the analysis is illustrated in Fig.|l|, 
where both the virtual and real photon processes are shown, as well as the region of track 
acceptance. 

2 Detector and Data Sample 

The data for this analysis were collected by the HI experiment in 1994, when HERA 
operated with 27.5 GeV positrons and 820 GeV protons, thus producing e + p interactions 
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Pt, max 



Figure 1: Illustration of diagrams for a DIS event and photoproduction event with high 
Pt particles in the final state. The sketch on the right side illustrates an inclusive charged 
particle pseudorapidity distribution in the 7p centre of mass system, with the range used 
to search for highpt tracks. For real photon interactions all contributing processes (direct, 
resolved, vector meson production, ...) are considered. 



at a centre of mass energy of y/s = 300 GeV. A detailed description of the HI apparatus 



can be found elsewhere |TI[ • The detector covers almost Att of the solid angle and consists, 
moving outwards from the interaction point, of tracking detectors, calorimeters and muon 
detectors. Along the beam line several other dedicated detectors, such as the luminosity 
system, have been installed. This analysis makes use of the luminosity system and the 
central tracking detector. 

The scattered positron was detected in the electron tagger of the HI luminosity system. 
The system consists of two TICl/TlBr crystal calorimeters, installed in the HERA tunnel, 
each having a resolution of a(E)/E = 0.1/y/E with E in GeV. The electron tagger is 
located at z = —33 m and the photon tagger at z — —103 m from the interaction point 
in the direction of outgoing positron beam. The energy scale of both calorimeters is 
known to better than 1.5% [I2|. The acceptance of the electron tagger is confined to the 



kinematical range of 0.2 < y < 0.8 and Q 2 < 0.01 GeV 2 . Here y is the fractional energy of 
the photon and can be well approximated by y = 1 — E' e /E e , with E e , E' e the energies of 
the initial and scattered positrons respectively. The photon tagger accepts photons with 
7i — 9 < 0.45 mrad. Here 9 is defined with respect to the incident proton direction. 

The charge and momentum of charged particles were measured by two coaxial cylin- 
drical drift chambers (central jet chamber, CJC), which cover the polar angular range 
15° < 9 < 165°. A superconducting solenoid provides a uniform magnetic field of 1.15 T 
parallel to the beam axis in the tracking region. The p t and 9 resolutions of the CJC are 
(Jpjpt ~ 0.009 • p t [GeV] © 0.015 and a e = 20 mrad respectively. 

Photoproduction events were triggered by a coincidence of an energy deposit larger 
than ~ 5 GeV in the electron tagger and one or more track candidates with p t > 0.4 



5 



GeV from the CJC trigger [13]. In order to reduce the non-ep background contamination 
and to ensure good reconstruction of the event kinematics, the following selection criteria 
have been applied for this analysis: 



• The event vertex reconstructed from the tracks of charged particles was required 
to lie within ±30 cm of the mean ^-position of the interaction point, corresponding 
to a 3a cut on the interaction region defined by the longitudinal size of the proton 
bunch. 

• In order to suppress random coincidences between the high rate of Bethe-Heitler 
events ep — > epj andp-gas background in the main HI detector, events were rejected 
if a photon with an energy E 1 > 2 GeV was detected in the photon tagger. This 
cut also reduces QED radiative corrections to the ep Born cross section 0. 

• The y range was limited to 0.25 < y < 0.70 to avoid regions of low electron tagger 
acceptance Jl]]. This corresponds to photon-proton collision energies of 150 < W = 
^ys < 250 GeV. 

• To ensure a good measurement of the scattered positron energy, a fiducial cut was 
applied to reject events with a positron detected close to the acceptance boundaries 
of the electron tagger. 



In total 1.8 • 10 6 events satisfy the above selection criteria, from a data sample which 
corresponds to an integrated luminosity of £ = 2.6 pb _1 . The residual background con- 
tamination was estimated to be (2.4 ± 0.5)%. 

The data were divided into several sub-samples according to the charged particle 
with the maximum transverse momentum, p^max; found in the pseudorapidity range 
1.1 < r]* < 3.1. Here pseudorapidity is defined in the photon-proton centre of mass 
system as 77* = — lntan(0*/2) with the polar angle 9* measured with respect to the 
photon direction. The pseudorapidity interval chosen corresponds to a region in the 
detector which guarantees a good track measurement in the CJC and a uniform track 
acceptance over the whole W range studied. Only those tracks were considered which 
have p t > 0.15 GeV, are fitted to the primary event vertex and further satisfy a number 
of quality criteria detailed in Ref. [H]. The systematic error due to the track selection was 



estimated by varying these quality criteria in a wide range and repeating the complete 
analysis. In total eleven bins in p^max have been defined, covering the Pt,max range from 
0.5 to 7 GeV. The bin sizes were determined to have adequate statistics and such that 
the size is always larger than four times the p t resolution in the bin to minimize migration 
effects. The values of the pt,max boundaries for the bins are indicated in Fig. 2. 

The track reconstruction efficiency, the contamination from secondary vertices and 
decays close to the primary vertex, and residual bin-to-bin migrations have been deter- 
mined from Monte Carlo studies. Events were generated with the PYTHIA model |15| 



and passed through the full HI detector simulation. All simulated events were subjected 
to the same reconstruction and selection procedure as the real data. The migration ef- 
fects were found to be small, between 10% and 20%, and the ratio of the number of 
reconstructed to the number of generated events to be close to one (±4%) in all p^max 
bins. 
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3 Results and Discussion 



In ep collisions, for any sub-class i of photoproduction events the energy dependence 
<7* (W) can be extracted from the corresponding differential ep cross section using the 



Weizsacker- Williams formula [15] for the photon flux F(y, Q ). After integrating over Q 



the corresponding relation reads: 



dN l 

ep -(l + 5 RC )£F(y) A(y) e\y) aUW), 



dy 



where N l ep denotes the observed number of events belonging to class i, L is the integrated 
luminosity, A{y) is the acceptance of the electron tagger and e l {y) is the efficiency of the 
trigger and selection criteria. The factor (1 + <5rc) is a correction for the cross section for 
QED radiation to obtain the ep Born cross section. 

The main systematic uncertainty, which so far prevents an accurate measurement 
of the energy dependence of the total photoproduction cross section at HERA, results 
from the large variation of the tagging acceptance A(y) (between 0.15 and 0.8 in the W 
range of interest), and which, in addition, strongly depends on the beam optics [Q, [T2 



However, the relative energy dependence for different photoproduction event classes can 
be determined with much better precision. Hence we consider the ratio 

[!J) o%(W) dNjdy e%y) ~ f&y) ' 1 J 

where e and e l denote the efficiencies for all events and for the events in class i respectively, 
and R is the uncorrected ratio. All efficiencies in (Q) cancel except for the trigger efficiency. 
Hence r\ is the ratio of the trigger efficiency for the sample i to that of the full sample. 

Applying this method the W dependence of the uncorrected ratio R for each bin i of 
Pt,max is shown in Fig. ^| together with a power-law fit of the form R oc W 2 ' x( - Pt ' m&x ' . 

In order to obtain the energy dependence of the partial cross sections, defined in a 
Pt,max bin with width 2 A, 



„ f'Pt. max+A 

a y P ( W ) = a i P ( W 'Pt,max) = / da Jp (W)/dp max ■ G^max 



in (0) the parameters X(pt,max) have been corrected for trigger efficiency r\(y) and account 
has been taken of the energy dependence of the total cross section in the denominator of 

® : 

) + Ajp (3) 

where X t (pt,ma,x) is the trigger efficiency correction and Ajp = 0.08 ± 0.03. The rather 
conservative error assigned to Xjp accounts for the invisible part of the total 7p cross 
section with the data selection of this analysis (mainly elastic and low-mass diffractive 
dissociation processes) The correction due to the trigger efficiency was obtained 
directly from the data, using a sub-sample taken with a trigger based solely on calorimetric 
quantities. Correlated systematic errors of the r*(y) largely cancel in the measurement 
of the slope parameter A 7P . The contributions of the systematic errors to the total error 
5A 7P , excluding the overall uncertainty in Ajp, are summarized in Table 1. The rightmost 
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Figure 2: The energy dependence of the ratio R in different bins of the Pt,max. of the 
charged particles produced in the pseudorapidity window 1.1 < r/* < 3.1. The curves 
represent a fit of the form R oc W 2 * A (Pt.max) motivated by the Regge formalism. For 
visibility the data of the three lowest pt,max bins are multiplied with the factors indicated 
in the figure. 



column gives the systematic error resulting from the precision of the W 2 reconstruction 
and the stability of the power-law fit to variations of the W 2 range used. 

The final corrected values of A 7P (p fimax ) are plotted in Fig.^|a against the scale which 
is chosen to be (2pi iinax ) 2 . The statistical errors dominate the precision of the \ lp (pt,max) 
measurement in the high-p t range (pt,max > 3 GeV). 

A fit to the data of the form A 7P (pt jmax ) = C ■ ln(4p 2 max ) + A with full errors except 
for 6\jp (full straight line on Fig.|3]a) gives a x 2 /ndf = 0.3 with C = 0.112 ± 0.007 and 
A = 0.031 ± 0.014. Note that when the scale 4p 2 max is multiplied with a constant, this 
does not change the slope C of the scale dependence, but only influences the value of A . 
Thus the change of A 7P (pt imax ) with the scale has physical significance. 

Fig.|^a shows also predictions of a longitudinal phase space model (LPS) and of the 
PYTHIA model ||15|| . The PYTHIA model is based on perturbative QCD for the hard 



scattering processes, but also includes a Regge inspired soft interaction component. For 
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Table 1: Values and the systematic error summary for the measurement of the slope 
parameter X lp {pt,ma,x) ■ The overall systematic error 5\p = 0.03 is not included in this 
table, but is included in Figfja. 



Bin 


rt,max 

[GeV 2 l 


\^r,(Vi mav) ± 5 X^ni StCtt) 

lP V-i E .lllcLA. / — ' V / 


Systematic errors 


track sel. 


trier f?er 


W 2 fit 


1 


0.36 


0.080 ± 0.006 


0.001 


0.02 


0.006 


2 


0.56 


0.121 ± 0.009 


0.001 


0.01 


0.008 


3 


0.80 


0.157 ± 0.007 


0.001 


0.01 


0.007 


4 


1.29 


0.211 ±0.008 


0.002 


0.01 


0.004 


5 


2.06 


0.267 ±0.012 


0.001 


0.01 


0.004 


6 


3.14 


0.334 ±0.015 


0.001 


0.01 


0.007 


7 


4.90 


0.344 ±0.022 


0.001 


0.01 


0.003 


8 


7.77 


0.407 ±0.031 


0.001 


0.01 


0.006 


9 


12.46 


0.484 ±0.051 


0.015 


0.02 


0.008 


10 


19.56 


0.565 ±0.079 


0.036 


0.03 


0.004 


11 


33.43 


0.603 ±0.105 


0.071 


0.07 


0.062 



the LPS model particles were generated with a uniform distribution in rapidity, taking into 
account energy-momentum constraints, and with a transverse momentum distribution 
da/dp 2 oc (p 2 + pl)~ n with p = 0.5 GeV and n = 2.9, tuned to describe the charged 
particle multiplicity and pt spectra of the data in the photon fragmentation region. As 
expected, phase space effects start to become important only at very large pt,m&x values 
and cannot explain the trend observed in the data. PYTHIA, on the other hand, describes 
the data well. In studying the different ingredients in this model it turns out that the 
main part of the rise can be explained by the integration of the leading order matrix 
elements over the available phase space. This is in agreement with expectations reported 
in |TJ]]. Since this integration is regulated by the quantity p 2 /s, with s the invariant mass 
squared of the hard scattering subsystem, the rise for a given pt,max value will depend on 
the W range studied. Parton showers, simulating next to leading order effects, as well as 
multiple parton-parton scattering, account for at most 20% of the measured A 7P (pt )max ). 
The effects of using different proton and photon structure functions are small. Hence the 
change of A 7P (pt ima x) is not associated with the growth of the gluon density in the proton 
at small-x, contrary to naive expectation. 

In Fig.|3|b the measurement of A in cr*°p oc (W 2s ) x dis(Q 2 ) j s shown for the virtual-photon 
proton cross section as function of the virtuality Q 2 [|J. For DIS the Q 2 dependence of 
Xdis(Q 2 ) is qualitatively described by perturbative QCD, as shown in Fig. 3b by the 
calculations for two quite different proton structure function parametrisations [18], 



Both are based on the DGLAP evolution equations [20[] but assume a different x-behaviour 
at the starting scale, Qq, of the evolution. At Q 2 ~ 4 GeV 2 the MRS-DO' parton densities 
are constant in x (A ~ 0), while for GRV they are singular at x = (A > 0). This 
difference is reflected in the predictions for the GRV and MRS-DO' parametrisations in 
Fig. 3b. However, both predict a significant rise of Xdis(Q 2 ) with Q 2 , due to the singular 
behaviour of the box diagram g — ► qq of the DGLAP evolution kernel In detail 

neither of the two calculations agrees very well with the data which are in between the 
predictions. Note that the MRS-DO' distributions were found to be in disagreement with 
the structure function data itself 0. Bl . 
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Figure 3: a) The scale dependence of the slope A 7P (pt imax ) in photoproduction (full 
points). The full line represents the result of a linear fit X-y P (pt,ma,x) = C ■ ln(4p 2 max ) + A 
to the data. The dashed line is the prediction from PYTHIA and the open points are 
the prediction from a longitudinal phase space model, b) The scale dependence of the 
slope \dis(Q 2 ) in deep inelastic data compared to different parametrisations of the proton 
structure function F 2 : full line GRV, dashed line MRS-DO'. For both figures inner error 
bars show the statistical errors; full error bars correspond to statistical and systematic 
errors added in quadrature. 



A qualitative similarity is seen between the W dependence of the cross sections for 
high-pt particle photoproduction and for deep inelastic scattering: they both show a 
significant rise of A 7 p(p timax ) and Xdis(Q 2 ) with increasing scale. The use of different 
scales in the study of this effect in the two processes prevents however a direct quantitative 
comparison. The range of the scale is varied from 1 GeV 2 to 100 GeV 2 and within 
that range the rise of the slope parameter is stronger for a yp (W, p 2 max )- In view of the 
discussion above the rise in both processes can be understood as follows: 

• Perturbative QCD : as argued above, the rise of the cross section for photoproduction 
events with a high p t particle can be explained in terms of parton-parton scattering. 
On the other hand, the rise of the slope parameter due to QCD evolution in DIS can 
be traced to the singular behaviour of the box diagram describing the elastic 7*p 
scattering amplitude. This diagram gives the contribution to the 7*p cross section 
of photon-parton scattering with gluon or quark radiation and can be expected 
to have similar properties as the diagrams for parton-parton scattering in high p t 
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photoproduction. It can be shown |22| that the photon-gluon fusion diagram leads 
to a rise of F 2 which is compatible with the HERA data . Hence, in view of Fig. 1, 
the cross section evolution with the scale (p 2 in photoproduction and Q 2 in DIS) 
can then be readily explained by leading order partonic scattering contributions. 

Regge approach: in Reggeon Field Theory, high energy particle interactions are 
described by pomeron exchange and A is related to the pomeron intercept. At high 
energies the unitarity corrections to the one-pomeron exchange, which effectively 
reduce the cross section value, become important. It has been shown |24| that these 
reduce the value of A. In both cases studied here the cross sections are far below the 
unitarity bound [53 and in a region where unitarity corrections become gradually 
smaller ||. Hence the cross section can rise faster with W for either increasing p\ 
or Q 2 . The data shows that this is indeed the case and the development of the rise 
is similar, albeit not identical, for both cases. 



In short, an interesting similarity is observed in the evolution of the energy dependence 
of a 7P (W, £>? )max ) and a^* p (W, Q 2 ), in qualitative agreement with leading order perturba- 
tive QCD and Regge Field Theory expectations. When decreasing the scale from large 
values towards 1 GeV 2 , which corresponds to a p^max of 0.5 GeV, no clear change is 
observed in the scale dependence of \ lp (pt,ma,x) ■ Hence the evolution towards the soft 
region appears to be rather smooth. 



4 Conclusion 



A significant change is measured in the energy dependence of the partial cross section 
a lp (W, p 2 max ) with increasing transverse momentum of the charged particles in the photon 
fragmentation region. This result cannot be described by a simple longitudinal phase space 
model. It is, however, well reproduced by leading order QCD and can be qualitatively 
understood in the framework of Reggeon Field Theory. A similarity is observed in the 
behaviour of the energy dependence of the <r lp (W, Pt imax ) an d cr*** (W, Q 2 ). This supports 
a common underlying picture for real and virtual photoproduction as well as a smooth 
and continuous transition from hard to soft processes. 
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